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Abstract. We theoretically and experimentally examine the effects of anharmonic 
terms in the trapping potential for linear chains of trapped ions. We concentrate 
on two different effects that become significant at different levels of anharmonicity. 
The first is a modification of the oscillation frequencies and amplitudes of the ions' 
normal modes of vibration for multi-ion crystals, resulting from each ion experiencing 
a different curvature in the potential. In the second effect, which occurs with increased 
anharmonicity or higher excitation amplitude, amplitude-dependent shifts of the 
normal-mode frequencies become important. We evaluate normal-mode frequency and 
amplitude shifts, and comment on the implications for quantum information processing 
and quantum state engineering. Since the ratio of the anharmonic to harmonic terms 
typically increases as the ion-electrode distance decreases, anharmonic effects will 
become more significant as ion trap sizes are reduced. To avoid unwanted problems, 
anharmonicities should therefore be taken into account at the design stage of trap 
development. 



I Author to whom any correspondence should be addressed. 



Normal modes of trapped ions in the presence of anharmonic trap potentials. 2 
1. Introduction 

Quantum control of the motional degrees of freedom of trapped atomic ions is among the 
most advanced in physics. It has an important role in realizing deterministic quantum 
information processing with trapped ions and offers a rich playground for the exploration 
of quantum-state engineering and quantum control. Examples of the latter include 
the production of Fock states, squeezed states, coherent states, and superpositions of 
coherent states analogous to the "Schrodinger's cat" thought experiment [U, El |3l IH |5] . 
The motional states of trapped ions have been a key ingredient for simulation of 
quantum-optical systems [6] and the Dirac equation [TJ [8] and for the realization of 
quantum walks in phase space [9l [TO] . 

For ion-trap quantum information processing pLlj, the collective motion of the ions is 
critical for performing deterministic multi-qubit quantum gates. Recent work includes 
the entanglement of up to 14 ions [12], and a two-qubit entangled state fidelity of 
99.3% [13]. In order to perform more complex operations, higher precision control 
of both the internal and motional states of the ions will be required, and scalable 
techniques will need to be implemented [H]. The high operation fidelities of around 
99.99% required to achieve fault-tolerance for quantum information processing [T5l 116] 
place stringent demands on the motional control. 

Most of the experimental work performed to date has used traps where the 
potential is harmonic to a very high degree. The resultant normal-mode frequencies 
and amplitudes have been extensively studied [171 [THl [19]. The relative importance of 
anharmonic terms depends primarily on their length scales relative to those relevant 
to the ions. In this paper, we consider two important effects that occur at different 
length scales. First, weak anharmonicity modifies the normal modes of multi-ion chains 
when anharmonic terms become significant over the length L of a chain of ions. In 
this case the ions' equilibrium positions are modified and the local curvature of the 
potential is different at the position of each ion. This modifies both the frequency and 
ion amplitudes of the normal modes of the chain, but the normal-mode eigenstates are 
still stationary states of the Hamiltonian. Second, strong anharmonicity plays a role 
when changes in curvature over the size of the motional wavefunction become significant. 
This creates motional frequency shifts as a function of the motional excitation and cross- 
coupling between different normal modes. The effects of strong anharmonicity have been 
observed when caused by anharmonicity in the Coulomb interaction between the ions 
[201 [21]; we extend these results to include anharmonicity in the trapping potential. The 
Coulomb anharmonicity plays a crucial role in the prediction of a temperature-driven 
structural phase transition in ion chains [221 123] , and anharmonic trap potentials will 
modify the effect and perhaps allow it to be tuned. Although we focus in this paper on 
anharmonicity as a perturbation to harmonic trap potentials, several experiments use 
strongly anharmonic potentials to split and recombine chains of ions [T¥l [2H [25l [26] and 
to engineer ion chains with nearly equal ion spacing [27] . 

Manipulating ions in small-scale traps that can be microfabricated is useful for 
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large-scale quantum information processing [HI [SSI [2S]; for achieving more intricate 
motional control [251 EOl EH [32l |33] ; and for coupling trapped ions to other quantum 
devices [HU [351 ESI EZ] • Anharmonic terms have the potential to scale unfavorably as trap 
sizes are reduced. Consider a trap with some characteristic length scale p. This length 
might be, for example, the distance from a trapped ion to the nearest point on the surface 
of a trap electrode. In an expansion of the trap potential VJ, the size of the nth term is 
proportional to d'^Vt/dz^ and thus it will tend to scale as p"""*"^ relative to the harmonic 
term. Moreover, some recent microfabricated traps employ a geometry in which the 
ion resides above a surface containing the electrodes [291 EH ESI HQ]. The intrinsic 
asymmetry in these surface-electrode traps can lead to large odd-order anharmonicities 
in the direction perpendicular to the electrode surface. While anharmonicities should 
tend to increase with decreasing trap size, the trap potential can still be engineered to 
suppress them, provided the trap is designed with sufficient degrees of freedom. 

To give a sense of the typical lengths involved, we note that a single ^Be"*" ion in 
the ground state of a 1 MHz harmonic potential has a root-mean-square wavefunction 
size of cr = 24 nm. Two ions in the same potential are separated by L = 9 pm. Recent 
trends in quantum information are towards increased numbers of ions to perform multi- 
qubit operations simultaneously on many qubits [121 HB 1121 SSI HI] , to accommodate 
sympathetic refrigerator ions along with those used to store qubits [HI [251 [SHI [27], and 
to encode logical qubits in multiple physical qubits [IHl HHl HT] . For chains containing 
eight ions, trap frequencies along the chain direction between 1 MHz and 5 MHz lead 
to a range of L between 36 /im and 12 /im. 

The paper is organized as follows. After a brief introduction to calculation methods 
for normal modes of trapped-ion crystals, we discuss weak anharmonicities, giving 
simple examples, experimental results, and methods for characterization. Next, we 
theoretically treat effects arising from stronger anharmonicities. We evaluate these 
effects with reference to a surface-electrode trap that has been used at NIST. In both 
sections, we discuss the implications of our results for precise control of trapped-ion 
chains for quantum information and state-engineering, focusing on multi-qubit quantum 
logic gates that use the ions' normal modes. Finally, we provide a short discussion of 
the susceptibility of anharmonic traps to electric field drift, after which we conclude. 

2. Calculating normal modes of trapped ion chains 

Normal-mode frequencies and amplitudes may be calculated from the classical 
Lagrangian equations of motion and quantized with the usual harmonic-oscillator 
algebraic formalism. We consider the case of ions each with charge g in a potential 
well. The energy of the ions is given as a sum of kinetic- and potential-energy terms 
E = T + U, where 



N 




i=l 
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Here mi, Zi denote the mass and position of the zth ion. The trap potential Vt includes 
both a mass-independent static potential and a mass-dependent pseudopotential arising 
from the radiofrequency confinement [2S]. In what follows, it is convenient to write the 
coordinates as 3A^ scalar parameters ^i, . . . , z^^ that produce the vectors Zi, . . . , z^v- 
The set of equations dU/dzi = give the set of equilibrium positions {z^} for the 
ions. For large this is done by numerical minimization, since the analytical expressions 
become complicated. In a Taylor expansion of the potential around these equilibrium 
positions, the leading term is at second order, which gives the symmetric Hessian matrix 

1 d^U 



H' = 
^■^ y/rftjnj dzidzj 



(2) 

As is typical [48j, we use mass- weighted coordinates; we indicate this transformation 
with a prime, = ^/rnlzi. This transformation allows us to write the kinetic energy T 
in a form that is independent of mass. 

The normal modes and their corresponding frequencies can be found by solving the 
Lagrangian equations of motion for the system. In this case, the relevant quantities are 
the displacements from equilibrium, Q = z[ — zf . Neglecting higher orders in U than 
those described by ([2]), the 3A^ equations of motion are 

ZN 

c; + 5^H^c; = o. (3) 

i=l 

Inserting a fiducial solution Q = Cj'e*^* gives a linear system of equations that can 
be diagonalized to yield the normal modes of the system, which are defined by the 
eigenvalues and eigenvectors of the matrix H'^j. The eigenvalues are equal to u"^, where 
Ua is the motional frequency of the normal mode a. The matrix of eigenvectors e'f allows 
us to express normal-mode coordinates as a function of the individual ion coordinates 
by use of 

i=l 

Since each normal mode acts as an independent oscillator, we can quantize them in the 
usual manner, writing the mass-weighted position operator as 

C = (^'a («a + «1) , (5) 

where a'^ = ^JhJi^Lj^ and a„ are the raising and lowering ladder operators. For 
the ith ion in mode a, the ground-state wavefunction root-mean-square size is 

= 4=er.i. (6) 



The quantized form for the ion's excursion from equilibrium is 



a=l 
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2.1. One- dimensional simplification 

For much of the paper, and in many experiments, we are primarily interested in the 
dynamics of motion along the ion chain. For an axial trap potential centred at z = 0, 
we can expand it as a power series: 



where A„ = {hu/ K'2Y^'^'^~^^ is a length used to parameterize the anharmonicity. We will 
consider potentials where the harmonic {n = 2) term dominates, as is typically the 
case in the experiments described below. In terms of the equation above, we assume 
I A/ An I <^ 1, where A is some relevant length scale such as the ion chain length L or 
the ground-state wavepacket size a. We assume no axial pseudopotential, so the series 
co-efficients do not depend on the ion mass. 

3. Anhctrmonic modifications to normal modes 

For anharmonic perturbations that appear on length scales comparable to that of the 
ion chain but still negligible on the scale of an ion wave-packet, the main effects are 
modifications of the normal-mode frequencies and amplitudes. For pairs of ions, these 
may be examined analytically, and we give examples of cubic and quartic potentials 
for pairs with equal and with unequal masses. The shifts seen are similar to those 
in longer chains, which are more easily analyzed numerically. We then examine these 
effects experimentally. We measure the anharmonic frequency shifts in chains of up to 
eight ions of equal mass. With mixed-species ion chains, we measure modified mode 
amplitudes and demonstrate a technique for nulling odd-order anharmonicities. 

3.1. Illustrative examples 

We start by giving two simple examples with equal-mass ions, which serve to illustrate 
the main effects that are observed in larger crystals, where numerical methods are more 
convenient. The parameter I — {q/87reoH2)^^^ is a characteristic length scale for the few- 
ion cases. For example, the distance between two ions in a harmonic well is L = 
The length of a chain of N ions in a harmonic potential scales as L ~ IN^-^"^ , obtained 
from a fit to numerically calculated chain lengths. In the equal-mass examples below, 
both cubic and quartic perturbations create a frequency shift that scales as (//A„)^. 
Wc consider the same perturbations for two unequal mass ions and show that now the 
frequency shift scales differently in the sizes of the perturbations, (//A„)"^^. In addition, 
the cubic perturbation depends on ion-order, while the quartic one does not, allowing 
the two types of anharmonicity to be differentiated. 

3.1.1. Example 1: cubic term with two ions of equal mass. Consider a potential of the 
form Vt{z) — {k,2z'^) (1 + z/Xs), and assume that |//'^3| ^ 1; that the cubic term is 
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small compared to the quadratic term over the length scale of the ion separation. The 
equilibrium positions z\ are found by solving the set of equations dU/dzi = 0, which 
has solutions (to second order in //A3) 

2" 



^± — 



"22/3 



It 



3 / 



25/3 A, 2V3 



(9) 



Expanding the potential about these equilibrium positions, we find that the 
eigenfrequencies are 



Us 



2qK2 


1 8 


m 




6qK2 


27/3 


m 



The corresponding eigenvectors are 

1 




3 j_ 

2573 A3 

^ ^ 25/3 A3 



and 





(10) 



(12) 



(13) 



In the limit where //A3 — ?■ these normal modes would correspond to the centre-of-mass 
and stretch modes. We see from these results that the corrections to the eigenvectors 
enter at lower order in //A3 than the corrections to the eigenfrequencies, as would be 
expected from general considerations of perturbation theory [59]. 



3.1.2. Example 2: quartic term with two ions of equal mass. The quartic term in the 
trapping potential has been studied for two ions in the context of separation of ions from 
a single- well potential into a double- well |SU], in which the quartic term was as strong 
as the harmonic term. Here we take an approach similar to that of the previous section 
and assume it to be small compared to the harmonic term. The potential is then of the 
form Vt{z) = /«2-2^ [1 + [z/XiY], with I//A4P ^ 1. In this case the centre of the ion pair 



remains at z 



0, and the equilibrium positions of the two ions are 

2 2^/3 / ' \ 41 
+ 



"22/3 



1 - 



3 • 21/3 



a; 



9 



(14) 



The normal mode frequencies are 
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Note that the frequency shift comes in as (Z/A„)^ for both the quartic and the cubic 
perturbations. For two ions, the normal-mode amplitudes remain the centre-of-mass 
and stretch modes, as is the case for a purely harmonic potential. For more than two 
ions, adjacent ions no longer see identical potentials, and the normal- mode amplitudes 
will shift from their harmonic values. 



3.1.3. Example 3: cubic term with two ions of unequal mass. We relax the equal mass 
assumption of the prior two cases, and re-examine the cubic term for ions with a mass 
ratio /X = mi/m2 (with zi < Z2). Because of our assumption of no axial pseudopotential, 
the only mass-dependence is in the kinetic energy, so the ion equilibrium positions are 
unchanged from ([9]). The normal- mode frequencies, however, are now 



2qK2 
nil 



1 + /i ± V/^ - l^+l 



1/2 



It 



1 



I 



28/3 



TA3 



.(17) 



Note that the frequency shift is now first-order in //A3 and depends on the ion order 
due to the first-order dependence on the sign of //A3. In the limit of a purely harmonic 
potential, (JT71) reproduces the results of (IHIET]. The normal-mode eigenvectors are 



v 



It 



1 



2^ ^/J^?^ l+'-i ^3 

3 1+M 4 / 

25/3 



-1 l+^'i -^3 



/ 



(18) 



Here, r± = [±(yU — 1) -|- ^^/JjF^^Ji^\^]/ y/JI. As for the equal-mass case, the higher- 
frequency mode (subscript -|-) involves the ions moving out of phase and the other 
mode involves in-phase motion. As before, the mode- amplitude corrections scale as 
//A3, though with mass-dependent coefficients. 



3.1.4. Example 4- quartic term with two ions of unequal mass. For unequal-mass ions 
in a potential with a quartic perturbation, the normal-mode frequencies are 

1/2 



X 



2qK2 
nil 

1 + 



^1 + /i ± -xZ/i^ - /i + 1 

1 t(1 + ^) + 7v^ 



t^ 



1 



3 ■ 24/3 



a; 



(19) 



V^/i^ - + 1 

In this case, the frequency shift is second-order in //A4, as for the equal mass ions. As 
one would expect from the symmetry of the potential, the frequencies do not depend on 
ion order. The eigenvectors are 



e_|_ ~ 



v 



Tl- 



1 ± 



2^ 7^:5^1+' 



A' 
A4 ) 

2 



(20) 



2^ ^-ri^^iT+^i 

Unlike the equal-mass case, the different masses break the symmetry, and there are 
quartic corrections to the mode amplitudes. 
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3.2. Experimental demonstrations 

Using various combinations of ^Be"*" and ^^Mg"*" , we demonstrate normal-mode frequency 
and amplitude shifts arising from weak anharmonic perturbations. We also present 
a technique that uses the ion order of an unequal-mass pair to suppress odd-order 
anharmonicit ies . 

3.2.1. Ion trap and potential wells. The linear radiofrequency Paul trap used here 
(described in ^25]) stores the ^Be"*" and ^^Mg"*" ions. The trap is formed from two 125 /^m 
thick wafers of alumina, with the electrodes made of sputtered gold and laser-cut 20 /zm 
vacuum gaps between adjacent electrodes. A top view of the top wafer in the relevant 
region of the trap is shown in figure [H Radial confinement of the ions is provided by a 
pseudopotential derived from a quadrupole potential {V^eak ~ 200 — 300 V) oscillating 
at a radio frequency (fij,^ = 27r x 150 MHz) and resulting in radial secular frequencies 
of approximately 12 MHz for a single ^Be"*" ion. The experimental region of the trap 
has five pairs of control electrodes that can be used to create axial potential wells (the 
trap has three additional pairs of control electrodes, but these are placed far from the 
experimental region and have only a small effect on the axial potentials in this region). 
The pairs of electrodes are arranged such that one of each pair is on the top wafer and 
the other opposes it on the bottom wafer [21]. The nominal voltage applied to the 
electrodes in one pair is the same; however a (typically small) differential component is 
also used to null out stray static electric fields at the pseudopotential zero as well as 
the effects of trap imperfections. In order to design potential wells for the ions located 
at a particular position, we simulate the influence of each electrode at that position 
using the Boundary Element Method (BEM) [52] . Since the fields from each electrode 
can be superposed, this allows us to design sets of voltages that produce harmonic or 
anharmonic potential wells at any point along the trap axis. It should be noted that 
while the simulation provides guidance, it relies on dimensions taken from photos of the 
actual trap, and thus its accuracy is limited. 

For the quantum information experiments performed in this trap [TU [251 126], we 
used two base sets of potentials. We periodically adjusted them with shim voltages 
or scaled them to adjust motional frequencies, but most experiments and waveforms 
began and ended with one of these sets. The first set produced a single potential well, 
positioned close to electrode 2, with negligible anharmonic terms along the axis of the 
trap. Despite our assumption above, slight electrode asymmetries can produce a small 
axial pseudopotential; therefore we designed this potential well to overlap its minimum 
position with the axial pseudopotential minimum, which we found experimentally. 
The voltage configuration used to do this was V*^^^ = {3.7,1.3,1.5,4.1,3.7} V, where 
the numbers indicate the voltage applied to electrode pairs 1 to 5. The simulated 
potential produced by these voltages is shown in figure [2]^a). The quantum information 
experiments required simultaneous trapping of ions in two experiment zones, for which 
we used the double- well potential shown in figure [2t^b). This potential used the 
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Figure 1. A top view of one of tiie gold-plated alumina wafers that combines with a 
second miderlying wafer to trap the ions. The picture encompasses the electrodes in 
the region used in the experiments. 

voltage set F^^) = {3.8,0.0,2.8,0.0,3.7} V, producing two potential wells 240 /xm 
apart. Both the single- well and the double- well potentials were designed to produce 
a minimum at the same position, -120 /xm in figure 121 with the same curvature, 
K2 = muj1/{2q) = 1.3 x 10^ V-m~^, corresponding to Uz = 27r(2.7 MHz) and / = 3.8 /xm 
for beryllium ions. From BEM simulations of the potential, we estimate the anharmonic 
terms of the double-well potential to be A3 = —200 fim and A4 = 250 /im. By 
comparison, simulations of the single-well potential give A3 = 10 mm and A4 = 1 mm. 




-200 -100 100 200 -200 -100 100 200 

position / |am position / i^m 



Figure 2. Simulations of the two axial trapping potentials. The origin of the position 
co-ordinate is in the center of electrode 3. (a) The harmonic-potential well, designed 
to have negligible anharmonic terms along the axis of the trap, (b) The double-well 
potential, which has significant anharmonicities in both wells. Zero corresponds to the 
centre of electrode 3. 

3.2.2. Control and readout of trapped ions. Control of the ions' internal states and 
motion is primarily through standard techniques described elsewhere [25| 135] . A 
11.964 mT magnetic field establishes a quantization axis aligned at 45 degrees to both 
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the trap axis and the vector normal to the alumina wafers. The ions' motion is initially 
Doppler laser cooled to temperatures low enough that the root-mean-square motional 
wavefunction size of each mode is less than 100 nm. For some experiments, we proceed 
with cooling one or more modes to near the ground state by use of Raman sideband 
transitions. The ^Be"*" ions' internal states are initialized by optical pumping to the 
Si/2 1-^ = 2,Mp = 2) hyperfine level, which we label as \i). To observe and Doppler 
cool the ions, we use a 313 nm cT"'"-polarized laser beam tuned close to resonance with 
the \i) ^ P3/2 \F' = 3, M'p = 3) closed cycling transition. State-dependent resonant 
fluorescence from the ions is detected on a photomultiplier tube (PMT) during a 200 /is 
detection window. For an ion initialized in \i), we observe an average of approximately 
eight counts on the PMT. 

Transitions among internal and motional states of the ^Be^ ions are induced by 
stimulated Raman interactions. These use two laser beams derived from the same laser; 
the pair of beams has a relative detuning equal to that of a hyperfine "carrier" transition 
or a hyperfine transition plus or minus a motional frequency ("sideband" transitions). 
The laser is detuned -70 GHz from the S1/2 ^ P1/2 transition. The geometry of the 
beams relative to the ions determines whether the lasers interact with only the internal 
states or with the motion as well. For motion-sensitive transitions, both beams are 
aligned at 45 degrees to the axis of the trap (one anti-parallel and one perpendicular 
to the magnetic field), with their difference vector Sk aligned along the axis, with a 
magnitude of \5k\ = 27r x \/2/X, where A = 313 nm. We typically utilize transitions 
between the states \l) = S1/2 \F = 2,Mf = 1) and |t) = S1/2 \F = 1,Mf = 0), whose 
frequency difference has no first-order sensitivity to magnetic field ^6]. In order to make 
use of this transition, we first transfer population from \i) — t- \D, again with a stimulated 
Raman interaction. Subsequent to driving Raman transitions on \D -H- |t), we transfer 
11) — 7- \i) and |t) — )• S1/2 1-^ = ^,Mp = —1), where the latter has no transitions that 
are close to resonance with the detection laser frequency and thus gives negligible 
fluorescence. 

To measure trap frequencies, we resonantly excite the motion of the ions. An 
oscillating "tickle" voltage is applied to one electrode of the trap, resulting in an 
oscillating electric field at the position of the ions. This excites motion when tuned 
close to resonance with a normal mode. Large motional excitations {n > 100) may 
be observed as a decrease in fluorescence of state \i), but motion-sensitive stimulated 
Raman transitions enable the detection of motional excitation corresponding to a single 
quantum. Each method will be described in detail below. 

3.3. Frequency shifts, homogeneous ion crystals 

One way to characterize the anharmonicities is to examine the axial frequencies of the 
in-phase motion of chains of ^Be"*" ions as a function of the number of ions in the chain. 
With the ions initialized in the internal state \i), we apply the tickle voltage for 500 /is 
and subsequently detect. A resonant force significantly excites the motion, resulting in 
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a drop in the ion fluorescence level (see for example [53]) as the ions' Doppler shifts 
approach the linewidth of the optical transition (for the fluorescence to drop to one 
half requires an excitation of approximately 1000 phonons for typical trap frequencies). 
For each frequency setting of the tickle voltage, we repeat the experiment 500 times, 
and record the average number of photon counts observed during the 200 fis detection 
period. We fit the results with a Lorentzian to obtain the normal-mode frequency. 

The resonant frequencies for the in-phase mode of motion of linear chains of one to 
eight ^Be"*" ions in the potentials created by V^^^ and V*^^^ are shown in figure [3l For the 
nearly harmonic well, the in-phase mode frequency changes by -0.4(1) kHz between one 
ion and eight ions, an effect that is negligible on the plotted scale. However in the case 
of the more anharmonic potential well, there is a frequency shift of —2.59(3) kHz/ion. 
Both the cubic and quartic terms produce a linear frequency shift per ion up to eight 
ions, so we are unable to distinguish these components in this type of experiment. 



2.66 

N 
I 

S 2.65 



2.64 




2 4 6 8 

number of ions 

Figure 3. The frequency of the centre-of-mass mode as a function of the number 
of ^Be+ ions. Triangles indicate the frequencies measured in the anharmonic trap 
given by potentials V^^^; circles indicate the frequencies in the harmonic trap given by 
potentials V^^\ The line is a fit to the anharmonic trap's frequency shift. Error bars 
are smaller than the symbols. 



3.4- Frequency shifts, inhomogeneous ion crystals 

To isolate only the odd-order components of the anharmonicity, we make use of 
an asymmetric ion chain, here consisting of one ^Be"*" and one ^"^Mg^ ion. We 
deterministically order the ions as ^Be^ - ^^Mg"*" or ^^Mg"*" - ^Be"^ (see below) . The length 
of the pair is L = 4.8 /xm. Odd-order anharmonicities produce different frequency shifts 
for different ion order, whereas even terms produce identical shifts. 

The mass difference between the ions and the mass dependence of the radial 
pseudopotential allow deterministic initialization of the ion order |5H |55l |56]. The 
re-ordering proceeds by first applying a radial electric field under continuous Doppler 
cooling so the ions remain close to their equilibrium positions. This field displaces the 
24]\/[g+ ion more than the ^Be"*" one due to the weaker pseudopotential experienced by the 
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heavier-mass ion. At a critical value of the electric field the axis of the two-ion crystal 
is normal to the trap axis. At this point, a differential shim voltage between upper 
and lower electrodes (e.g. electrodes 2-upper and 4-lower raised relative to 2-lower and 
4-upper) can twist the axis of the static potential, breaking its symmetry relative to the 
orientation of the two ions, and pushing the ^^Mg"*" to a different value of z compared to 
the ^Be"*" ion. By subsequently reducing the field and removing the twist, the ions are 
left in the desired order. 

After ordering the ions, the frequency of the in-phase normal mode of the two ions 
is measured by use of the tickle method. In the case of two ions of unequal mass, 
both axial modes of motion involve modulation in the distance between the ions. This 
means that anharmonicity in the Coulomb interaction can produce frequency shifts on 
both modes as a function of motional excitation [20] . These shift the deduced motional 
frequency, and can prevent a single frequency drive from exciting the ions to energies 
high enough to reduce fiuorescence. 

To reduce the excitation energy at which a signal can be observed, we can observe 
weak motional excitation using Raman transitions. We first initialize both axial modes 
of the ion chain close to the ground state of motion by use of Raman sideband cooling 
[57] . After applying the tickle, we probe the motional excitation of the ^Be"*" ion by 
resonantly driving the IJ,) ^ |t) carrier transition with the motion-sensitive Raman 
beams. We choose the drive duration such that an ion in the ground state would make 
a full transition from from \\) — )■ |'|'), which would leave a subsequent detection "dark". 
If the motion of a normal mode is excited to state |n), the transition rate is reduced 
by a factor given by the matrix element (n| exp(z?7(a + a))) |n), where rj = Ska is the 
Lamb-Dicke parameter [35] . This results in incomplete population transfer, resulting in 
the detection of fiuorescence in the subsequent detection window. For the in-phase mode 
of motion of a ^Be"*" -^^Mg"*" pair, rj = 0.18, and the Rabi frequency of the transition 
is reduced to approximately half the ground-state rate for n = 17. A more precise 
description involves averaging the transferred population over the Fock state distribution 
of a coherent state, however the current estimate is sufficient for our present purposes. 

Data from frequency scans of the tickle voltage applied near the in-phase axial mode 
frequency for two configurations of ions in the anharmonic potential well characterized 
by V^'^^ are shown in figure HI The data sets are each fitted with a Lorentzian, and the 
frequency of the normal mode extracted. We measure a difference of 20.8(2) kHz between 
the mode frequencies for the two configurations. If the only odd-order component is a 
cubic term, this corresponds to A3 = —230 fim. By comparison, a similar experiment 
performed in the more harmonic well characterized by V^^^ gives a frequency shift of 
less than 0.03(2) kHz, corresponding to IA3I > 150 mm. 

An even more sensitive method of observing motional excitation is to probe on a 
motional sideband of the |t) transition. If the motional mode is in the ground 

state the motion-subtracting sideband cannot be driven. However, when the motion is 
excited, population transfer is allowed between \l) and |t)- To maximise the population 
in n = 1 (for which the Rabi frequency is set to make a full transfer), we would require 
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f/MHz f/MHz 

Figure 4. Frequency of the in-phase mode for ions in the order (a) ^^Mg"*" -^Be"*" and 
(b) ^Be"*" -^^Mg^ in the anharmonic weU situated at -120 /im in figure O 

the tickle to excite the motion to a coherent state j58j with a 1. A coherent state of 
this size would result in half the population being transferred from \\) to W). We use 
this method for nulling anharmonicity as described in section 13.51 

3.4-1 ■ Measurement of modified ion amplitudes. When driving mot ion- sensitive Raman 
transitions, the amplitude of motion of the ion in the normal mode affects the modulation 
index of the light, thus affecting the Rabi frequency at which the mode is driven. 
Mathematically, for ion j, the transition matrix element for the resonant transition 
between n) and n + 1) is proportional to (n| exp[i(5/c aj{a+a))] \n + 1) [35], where 
Gj is given in ([6]) and is proportional to the ion's motional amplitude through the mode's 
eigenvector. For aj <^ 1/Sk, this expression reduces to 6k Oj\/n + 1, showing direct 
proportionality to the amplitude. 

In the experiment, we use a four- ion chain containing two ^Be"*" and two ^^Mg"^ ions, 
which are initialized in the order ^Be"*" - ^^Mg"*" - ^"^Mg"*" - ^Be"^ prior to each run of the 
experimental sequence [25]. Only the beryllium ions interact with the 313 nm Raman 
light fields. The length of the chain is L = 10.8 /im. 

To equalize the intensity of each Raman light field on the two ^Be"*" ions, we 
separately observe the AC Stark shift of each beam. Either beam alone cannot drive a 
transition between and but the light causes a phase shift that we monitor with 
a Ramsey-type interference experiment [59] where one beam is applied between the two 
Ramsey pulses. By equalizing the rates at which these phases evolve for both ions, we 
equalize the relative electric field strength of each light field at each ion to better than 
2%. 

Once the light fields are equalized on the ions, we drive motional sidebands 
of the four-ion chain. The two highest-frequency axial modes have frequencies of 
/s = 5.5 MHz and /4 = 5.7 MHz and are particularly sensitive to the cubic term in the 
axial potential. In the nearly harmonic potential they have normal-mode eigenvectors 
e'3 = (0.629,-0.322,-0.322,0.629) and = (0.532,-0.465,0.465,-0.532). In the 
more anharmonic potential, the amplitudes are significantly different. This becomes 




Figure 5. Experimental results and numerical simulations of sideband flopping on the 
two Beryllium ions for (a) the third and (b) the fourth axial mode of a 4-ion ^Be^ - 
^*Mg^ - ^''Mg"^ -^Be"*" chain in the anharmonic trap. The amplitude is related to the 
population of ions in the |t) state, A = P(tt) + [^'(ti) + ^(it)]/2, and is proportional 
to the total fluorescence from both ions. 



obvious when the internal state populations are measured as a function of sideband 
drive duration, where beating behavior is observed due to the different Rabi rates 
of each ion. Results of such measurements are shown in figure [51 Also shown are 
simulations of sideband drives for two ^Be"*" ions with a ratio of amplitudes of ion motion 
{Ra = Wai/^'ai\) ^3 = 0.625 and = 0.500 (these values were chosen by eye to 
best fit the data). The simulation curves are obtained by integrating the Schrodinger 
equation, and subsequently adding a phenomenological decay to the Rabi oscillations 
in order to account for both motional and internal-state decoherence. Both simulations 
use the same value of the carrier Rabi frequency. For the cubic term of A3 = —230 fim 
obtained from the frequency shifts of the ^Be^-^^Mg"*" in-phase mode, we would expect 
these ratios to be -R3 = 0.644 and R4 = 0.499. The eigenvectors in this case are 
e'3 = (0.474,-0.167,-0.452,0.736) and = (0.686,-0.531,0.359,-0.342). Though 
the agreement between the two values of R^ seems good, choosing R3 = 0.644 produces 
a noticeable mismatch between data and theory. The reason for this discrepancy is 
currently not understood. 

Unequal amplitudes of motion are generally undesirable for performing multi-qubit 
logic gates on trapped ions. For the gate methods in common use today [HI |13l [6OI [61] , 
driven ion motion during the gate adds a phase conditioned on the qubits' joint state. In 
these gates, the unequal mode amplitudes can lead to a longer required gate duration, 
thereby increasing the probability of error from off-resonant photon scattering [62] . 
Additionally, this phase could be different for \'fD and due to the different ion- 
amplitudes, which increases the complexity of gate calibration. 

3.5. Nulling odd-order anharmonicities 

For the purposes of optimizing quantum logic gates on multi-ion crystals, it is desirable 
to be able to tune out the anharmonicities. In the trap used here, we cannot do this 
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Figure 6. lon-order-dependent frequency shift for the potential weh cfose to electrode 

2 in the double-well potential. The (a) in-phase and (b) out-of-phase mode frequency 
between the ^''Mg"'" -^Be^ and ^Be"*" -^^Mg"*" configurations as a function of electrode 

3 voltage, which is plotted as an offset from the simulated value. 

while maintaining a double-well potential if the positions of the minima of the two 
potential wells are fixed to be the same as the minima of the set of potentials V^^^ 
(the number of independent electrode voltages is insufficient to satisfy all the required 
constraints). However, if we remove this position constraint, it is possible to produce 
two wells with zero odd-order anharmonicity (though it seems that one voltage should 
not be sufficient to meet the constraints of the two wells, in practice the symmetry of 
the electrodes with respect to the well positions is sufficient to meet both constraints). 
We perform this optimization by observing the ion-order dependence of the in-phase 
and out-of-phase mode frequencies of a ^Be"*" - ^^Mg"*" pair as a function of the voltage 
V3 applied to the control electrode positioned between the two potential wells. The 
potentials used here are slightly modified from V*^^^ and are optimized in simulation 
for negligible cubic anharmonicity in both wells. The nulling procedure allows in-situ 
adjustment to account for differences between simulation and experiment. 

Data from such an optimization are shown in figure El which shows frequency 
differences when reversing ion order for the potential well close to electrode 2 in the 
double-well potential. The data show that the value of V3 where the in-phase mode 
shows zero frequency shift on inverting the ion order is different from the value of 
V3 where the same condition is met for the out-of-phase mode. This discrepancy is 
consistent with a small pseudopotential gradient along the axis of the trap, indicating 
that the adjusted position is no longer at the axial pseudopotential minimum. Since the 
pseudopotential has a different strength for both ions, a gradient pointing from ^^Mg^ to 
^Be^ will result in the spacing of the ions being reduced (and vice versa). This tends to 
increase the energy of the out-of-phase mode, but has a smaller impact on the in-phase 
mode. From the shift of the out-of-phase mode at the point where the in-phase shift is 
nulled, we can calculate a value of the pseudopotential gradient of 0.2 eV-m~^. 

To choose which value of V3 to use to minimize the ^Be^ - ^^Mg"^ - ^^Mg^ - 
^Be"^ sideband amplitude imbalance described in the previous section, we numerically 
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calculated the amplitude imbalance with the pseudopotential gradient fixed. We find 
that the parameters for which the amplitudes are balanced are close to those for which 
the ^Be"*" - ^^Mg"*" in-phase mode shift is minimized. For this condition, the resulting 
amplitudes for the four ion chain are calculated to be e'g = (0.631, —0.323, —0.322, 0.628) 
and = (0.530, -0.465, 0.467, -0.533). 



4. Strong anharmonic effects — modifications to the normal-mode picture 

As anharmonic effects increase to the level where individual ions sample the 
anharmonicity during excursions about equilibrium (that is cr/A„ is no longer negligible), 
the independent normal-mode picture for ion motion is no longer valid. The 
anharmonicities lead to coupling between the normal modes. The principal effects of 
this coupling are mode frequency shifts as a function of the number of quanta in both 
the motional mode of interest and other motional modes. Anharmonic effects of this 
kind were observed in [20] where the mode cross-coupling was due to anharmonicity in 
the Coulomb interaction between two ions. A variety of couplings due to the Coulomb 
anharmonicity are also discussed in [63]. In what follows, we extend these results by 
considering shifts arising from an arbitrary applied trapping potential. As a concrete 
example, we evaluate this in the context of a surface-electrode trap with the ion trapped 
30 /zm above the electrode plane. 



4.1. Perturbation theory 

If the harmonic term still dominates the potential energy, we may treat anharmonicity as 
a perturbation. Our goal is to calculate the frequency shift A/^ of a given mode (labeled 
Z in this section) for the transition between nz and nz + ^- The shift may depend on 
the occupation of that mode uz as well as that of all other modes {ua}. Because of 
the potential dependence on other modes, we again consider all three dimensions, as 
we did in section O and revert to the notation used in that section. The cubic and 
quartic terms in a Taylor expansion of the potential energy U will involve sums over the 
coefficients 

1 1 d'^U 



'■^'^ 3! ^minijnik dzidzjdzk 
1 1 d'^U 



.'(4) 



(21) 
(22) 



4! ^mimjmkmi dzidzjdzkdzi 

where the indices refer to ions not normal modes. Both the Coulomb interaction and 
any anharmonic terms in the trapping potential contribute to these coefficients. 

Because we are interested in perturbations to the normal modes of the system, we 
transform these tensors into the normal-mode basis in a manner analogous to the vector 
case of ([3]): 

Gfl =<^;.<E^'^'<4? (23) 
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i,j,k,l 

The inclusion of the a' coefficients is for notational convenience. In this notation, the 
cubic and quartic terms of the Taylor expansion are 

U'^'^ = E ^SU^- + + + (25) 



^(4) 



a,/3,7,<5 



(26) 



Here, we have written the position operator ('^ in terms of raising and lowering operators 
as in (111). 

In ffist-order perturbation theory, all odd-order contributions vanish exactly and 
the leading contribution is from the quartic term: 

AEi {na,nis,n^,ns) = {na,np,n^,ns \ U^^^\na,np,n^,ns) . (27) 

Only terms with coefficients G'^^zz ^^"^ ^'alzz have nonzero contributions to the 
uz nz + 1 frequency shift. 

The cubic terms give the leading contribution in second-order perturbation theory, 

|2 



AE2 {na,n,3,n^ 



\{na,nfi,n^ |[/(^) | n^, n/j, n^) | 



^ h [Uja {Ua - Ua) + W/3 (n^ - Up) + UJ., {u^ - fl^)] ' 



(28) 



These terms can potentially lead to a number of resonances. For example, it is possible 
for two modes to have a resonance allowing the destruction of two phonons from one of 

'(3) 

the modes and the creation of one phonon in the other. The G^^^, term in the potential 
creates such an interaction. In this discussion, we assume we are detuned from all such 
resonances so that perturbation theory remains valid. As ion numbers increase, mode 
density will increase and it will be harder to avoid these resonances. 

After evaluating the ladder-operator algebra for the leading perturbation-theory 
contributions above, we find that the frequency shift for the uz nz + 1 transition is 

/iA/z({n„},nz) = AEi({n«},nz + 1) - /\Ei{{n^},nz) 

+ A^2({n„}, uz + l)- AE2({n,}, uz) 
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Figure 7. Diagram of (a) a microfabricated surface-electrode trap and (b) a closer 
view at the trapping region. Colors differentiate the electrodes, which are marked with 
their bias voltages. The radiofrequency electrodes are at DC ground. A single ion 
would be trapped near the circle in part (b) and 30 /xm above the surface; two ions 
would lie approximately L — 4.5 /im apart with their axis parallel to the electrode 
axes. The intrinsic asymmetry of the electrodes relative to an ion's position leads to 
anharmonic perturbations in the trapping potential. 



72 



3Af 3Af 
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(29) 



4-2. Example case: a microfabricated surface trap 

One of the smaller traps that has been used at NIST has the electrodes lying in a plane 
and ^^Mg^ ions trapped 30 fim above that surface (a nearly identical trap was used in 
|64]). Figure [7] shows the trap geometry. The trap is operated with single- ion secular 
frequencies for the radial modes of /i = 7 MHz and /2 = 5 MHz, and an axial frequency 
of /a = 1.8 MHz. The asymmetry of the trap electrodes with respect to the position 
of the ion create cubic perturbations in directions out of the plane of the trap that 
dominate the third-order tensor G'^^\ 

Using fl29|) . we calculate a cross- coupling matrix x, which relates mode frequency 
shifts to excitation through 



A/, 



3N 



(30) 
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For a single ion in this trap, 
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(31) 



The largest shifts are a few parts in 10'' per quantum. 
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Multi-qubit gates in quantum information processing typically involve multiple ions 
in the same trap zone. The couplings among the six normal modes of two trapped ions 
for the trap of figure [7] are 
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The tensor is ordered such that the highest-frequency normal mode is at the top (index 
1), and the lowest is 6. Modes 5 and 6 are axial, and modes 1-4 are radial. The modes 
where the ions' amplitudes have opposite sign (the out-of-phase modes) are 2, 4 and 
5, while the in-phase modes are 1, 3 and 6. For comparison, if only the contribution 
from the Coulomb interaction is included (that is, if the trap potential were perfectly 
harmonic, having the same single ion trapping frequencies as above), the shifts are 
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On comparing these matrices we observe that in this trap the Coulomb anharmonicity 
still creates the largest mode-couplings, although some terms arising from the trap 
potential approach the strength of the Coulomb contribution. However whereas the 
Coulomb anharmonicity principally affects only modes of motion where the oscillation 
of the two ions has a sizeable differential component to the motion, the trap potential 
anharmonicities couple all modes. 

4.3. Coherence during quantum state manipulations 

The frequency shifts described in the previous section can become problematic both for 
motional state engineering [H [21 |3l HI [TOl [25] and in quantum information processing, 
where the motion of the ions is used in multi-qubit gates. To give an idea of the 
level at which cross- coupling from anharmonic terms could impede control in typical 
experiments, we provide two examples: the coherence of a superposition of motional 
Fock states and loss of fidelity in a two-qubit gate. 

Consider the superposition {\Oz) + \nz))/ for a mode Z, and assume for simplicity 
that all other "spectator" modes can be described by a single Doppler-limit temperature 
Tjj [65]. The coherence between the two motional states is given by the off-diagonal 
element of the density matrix, which evolves in time because of the anharmonic 
frequency differences. Taking the thermal average of the motional excitation of spectator 



Normal modes of trapped ions in the presence of anharmonic trap potentials. 



20 



modes yields 
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For example, if the superposition is between the ground and first excited states of 
the highest-frequency radial mode (mode 1) of a single ^^Mg'''ion in the surface trap 
presented in the previous section, and the other two modes are cooled to To = 0.7 mK, 
the coherence decays to 1/2 in 40 ms. Because the cross-couplings are themselves 
coherent, C recovers to approximately 80 % at l/xi2 ~ 0.4 s and in principle almost fully 
recovers at l/xis ~ 25 s. Such recoveries have been observed in other experiments pU] . 
For a motional superposition with (|0) + |10))/-\/2, the decay is ten-times faster. 
Recent experiments in other traps have produced superpositions of motional states 
with occupations greater than n = 100, which were manipulated over time scales of 
milliseconds [10]. If similar experiments were performed in the surface trap described 
here, the influence of the mode cross-coupling would need to be taken into consideration. 

All deterministic multi-qubit gates in trapped-ion quantum information processing 
rely on transient excitation of motion conditioned on the ions' internal (qubit) 
states [601 [nil [6ll [66l [67] . At present, typical implementations involve cooling the ions 
to the ground state of one or more modes of motion prior to implementing the gate. 
To increase the processing speed, it would be advantageous to perform high-fidelity 
quantum logic gates after only Doppler cooling. This is possible when the ions are well 
within the Lamb-Dicke regime and such gates have been demonstrated using lasers [68] 
and magnetic- field gradients [61]. In these cases, all the motional modes are thermally 
occupied with non-zero n, and anharmonic effects need to be considered. The gate 
methods used in these demonstrations are different versions of a common type, which 
make use of transient motional excitation by an internal-state-dependent force. This 
force may be applied optically [601 ED |69] or via microwaves [HH [70] , where the form of 
the state-dependence can be tailored by control of the driving fields. 

Here for simplicity, we consider the approach taken in [61], which assumed a 
Hamiltonian 

H = hQ{t)S, cos {ujt)r]{ae-'''^' + a^e''^^*) (35) 

where fl(t) is related to the laser fields used to produce the gate, w is a drive frequency 
typically near uz, is the sum of the Pauli Z operators |t) (tl ~ li) (il acting on each 
ion's spin, t] = \6kai\ = \6ka2\ is the Lamb-Dicke parameter for the interaction between 
the light fields and the two ions involved in the gate, and and a are the creation 
and annihilation operators which act on the mode chosen for the gate. For simplicity, 
in what follows we assume that fl{t) = Q for < t < tg and is zero otherwise. If 
Icj — w^l ^ I I, I a; I we can make a rotating wave approximation with respect to 
the motional frequencies, resulting in an evolution operator for the system given by 

U{t) = D{a{t)S,) exp{i^{t)S'^^) (36) 
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where -D(/3) is the motional state displacement operator exp(/3d''' — [3* a) ^58j, and 

a{t) = -%e-"^''^sm{6tl2), (37) 



m = ^ [sin(5t) - 6t] , (38) 
with 6 = u — ujz- 

A common method for characterizing the performance of such a gate is to examine 
the fidelity with which the entangled state iV'ideai) = (|i4) — ^ |tt)) / "\/2 is produced from 
the state |||) when the force pulse is applied in the first half of a spin-echo sequence 
and f2 = (5 and tq = 2tt/6. Ideally, these values result in a^Tc) = and ^{tg) = vr/2. 
For general Q;(r), $(r), the fidehty for producing I'j/'ideai) at the end of the spin-echo 
sequence is given by [H] 



F = KV^ideailplV^ideai)! = ^ + ^e-^'^W' + leH^WI^/^ siu [$ (t)] . (39) 

o o Z 

In the presence of anharmonicity, errors can enter into this gate in three ways. 
First, if the motion begins in a distribution of states, each state will have a different 
mode frequency and therefore a different detuning. Thus, in general the motional state 
will not return to its initial position at the end of the drive. For a detuning 6{l + e) with 
fractional error e ^ 1, this results in a{TG) — 7re. The resulting residual entanglement 
between the internal states and the motion reduces fidelity when the motional degree of 
freedom is traced out. Second, the different detunings cause each initial motional state 
to enclose a different area of phase space and thus to acquire a different phase. Again 
using a detuning error parameter e, we find that ^{tg) — (7r/2)(l — 2e). These phases 
must be averaged over using the distribution of initial motional states. Finally, the 
transient motional excitation itself will give a time-dependence to the detuning as the 
excitation increases then decreases during the gate. For the small excitations typical for 
present multi-qubit operations, this third effect is significantly smaller than the other 
two, and we ignore it in below. 

Let us assume that the gate is performed with the duration chosen such that 
tg = 2it/6 for the detuning corresponding to all modes in the ground state. For ions 
which start in an incoherent motional state distribution with mean quantum numbers 
in the set {n^}, the above expression then leads to a fidelity 

^ = 1 - ^ XI XzaXzpfic^ni^ + xlaK J • (40) 

The bars indicate an average over the incoherent distribution of n. For a thermal 
distribution, the mean-square n in the last term can be rewritten as = na{2na + 1). 
Failure of the motional states to return to their original position contributes two-thirds 
of the infidelity, with the other third due to the distribution of accumulated phases. 

Consider two ions in the surface trap described in the previous section, and a two- 
qubit gate that makes use of the radial rocking mode with frequency /2 = 6.8 MHz 
(as used in [64j). Implementing such a gate with all modes cooled to the Doppler-limit 
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temperature of = 0.7 mK produces an infidelity of 1 — F = 4 x 10~^[27r(l kHz)/5]^ 
if all other aspects of the gate are perfect. For multi-qubit gates performed in a 
similar manner on larger numbers of ions [I2l |l2], the increased number of motional 
modes will reduce the fidelity even further. Thus for detunings 5 > 27r(20 kHz), the 
infidelity is below the level of 10~^, which is often estimated as a requirement for fault- 
tolerant quantum information processing. Since the cubic terms in the potential scale 
unfavorably with reduced trap size, anharmonicity may become a significant source of 
error as traps become smaller. Various approaches may be used to implement a gate 
in a manner that suppresses these errors. For example, though the detuning in the 
above example is 5 = fi, by use of two pulses — one performed in each half of the spin- 
echo sequence [71] — the detuning could be increased to 5 = y/2^l, thus reducing the 
error. This approach also has the advantage that the residual displacement arising from 
errors in detuning can be arranged to be of opposite sign for the two pulses, and hence 
is canceled out. This is a special case of the more general composite-pulse schemes 
recently suggested by Hayes et al. [72], that has been used in gate operations in several 
experiments [HI [251 12H1 SSI [71] . Pulses with non-square shapes may also help in this 
regard [75]. This method can be extended by increasing the detuning and executing 
more, but smaller, loops in phase space [H]. At ion-electrode distances of approximately 
30 /im, the coherence of motional-state superpositions is currently limited by anomalous 
heating [TSj [761 [^ ? which is larger than the anharmonic infidelity. 

5. Sensitivity to electric fields 

The range of effects described above focus on problems related to anharmonic trapping 
potentials. In addition to these, a practical consideration is that anharmonicities 
introduce a dependence of the secular frequencies on the position of the ion. This means 
that uniform electric fields can displace the ions and cause the trap frequencies to shift. 
For the two-layer trap described in section 13.21 and a cubic term with A3 = —230 fim 
(as observed in section 13. 3p . a field of 2 V-m~^ causes a fractional frequency shift of 
~ 10~^ (corresponding to 3 kHz in that example). For multiple experiments [T^ l25l [26] 
performed using this trap, we observed that the trap frequency stability of the harmonic 
trapping potential was better than that of the anharmonic potential. This suggests that 
frequency instability may have been caused by slow fluctuations in stray electric fields, 
but this correlation has not been characterized. 

6. Conclusion. 

Anharmonic trapping potentials give rise to a number of effects that should be taken 
into account when performing quantum state engineering tasks in ion traps. As trap 
sizes are reduced, these effects become more significant. By building traps with intrinsic 
geometric symmetry, it should be possible to minimize odd-order anharmonicities. This 
consideration should be made at the trap design stage. In practice, the fabricated 
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trap is unlikely to be an exact realization of any design. For this purpose, it is 
desirable to include a sufficient number of independently controllable electrodes to allow 
anharmonicities to be canceled in situ. In this case, methods such as those given in 
section 13.51 can be used to null out unwanted terms. In addition to undesirable effects 
due to strong anharmonic trap potentials, the advent of small traps where these terms 
can be engineered might also enable novel state-preparation schemes |63] and tunable 
phase transitions [22l [23] . 
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